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Abstract 

Classical persistent homology is not tailored to study the action of transformation groups that are different 
from the group Homeo{X) of all self-homeomorphisms of a topological space X. In order to obtain better 
lower bounds for the natural pseudo-distance da associated with a group G C Homeo{X), we need to 
adapt persistent homology and consider G-invariant persistent homology. Roughly speaking, the main idea 
consists in defining persistent homology by means of a set of chains that is invariant under the action of 
G. In this paper we formalize this idea, and prove the stability of the persistent Betti number functions 
in G-invariant persistent homology with respect to the natural pseudo-distance da- We also show how 
G-invariant persistent homology could be used in applications concerning shape comparison. 

Keywords: Natural pseudo-distance, filtering function, group action, lower bound, stability, shape 
comparison 

2010 MSG: Primary 55N35, Secondary 68U05 



1. Introduction 

In many applicative problems we are interested in comparing two M^^-valued functions defined on a 
topological space, up to a certain group of tranformations. As an example, wc can think of the case of 
taking pictures of two objects A and B from every possible direction and comparing the sets of images we 
get. In such a case each image can be seen as a point in M*^, and our global measurement as a function 
(f : S"^ ^ M.^ , taking each direction (represented by a point in S"^ C K"^) to the picture wc get from that 
direction. In this case the position of the examined objects cannot be predetermined but we can control the 
direction of the camera that takes the pictures. As a consequence, two different sets of pictures (described by 
two different functions (p,tp : S'^ ^ M'^) can be considered similar if an orientation-preserving rigid motion 
g of exists, such that the picture of A taken from the direction of the unit vector v is similar to the 
picture of B taken from the direction of the unit vector g{v), for every v G S^. Formally speaking, the two 
different sets of pictures can be considered similar if infg£j^(52) max^£52 \\'fi{v) — i/^idiv))]]^ is small, where 
R{S^) denotes the group of orientation-preserving isometrics of S"^. 

The previous example illustrates the use of the following definition, where G°(X,]R'^) represents the set 
of all continuous functions from X to R*^. These functions are called k- dimensional filtering functions on 
the topological space X. 

Definition 1.1. Let X be a triangulable space. Let G be a subgroup of the group Homeo{X) of all 
homeomorphisms f : X ^ X. The pseudo-distance da ■ C"{X,R'') x C°{X,R'') R defined by setting 

dG{<P,''P) = inf max||v?(a;)- V(5(a^))lloo 

is called the natural pseudo-distance associated with the group G. 
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The previous definition generalizes the concept of natural pseudo-distance studied in [Ol [HI [3 HI [10] to 
the case G ^ Homeo{X), and is a particular case of the general setting described in [TT. The case that 
G is a proper subgroup of Homeo{X) is also examined in |2j E], and in |12j for the case of the group of 
diffeomorphisms (in an infinite dimensional setting). 

The pseudo-distance da is difficult to compute. Fortunately, if G = Homeo{X), Persistent Homology 
can be used to obtain lower bounds for (Iq. For example, if we denote by Dmatch the matching distance 
between the n-th persistent Betti number functions pf^ and of the functions (p and tpi we have that 

{'P,ip) (cf- [HIS])- For more details about Persistent Homology we refer the 

reader to 011]. 

A natural question arises: How could we obtain a lower bound for (Iq in the general case G 7^ Homeo{X)7 
Does an analogue of the concept of persistent Betti number function exist, suitable for getting a lower bound 
for dc? Since dHomeo(x) V') ^ dciip, i^), one could think of using the classical lower bounds for the natural 
pseudo-distance dHomeo(x) in order to get lower bounds for the pseudo-distance da- Before proceeding we 
illustrate an example, showing that in many cases this choice is not useful. 

Example 1.2. Let us consider an experimental setting where a robot is in the middle of a room, measuring 
its distance from the surrounding walls by a sensor, for each direction. This measurement can be formalized 
by a function ^ : 5^ — > M, where ^{v) equals minus the distance from the wall in the direction and verse 
represented by the unit vector for each w g S"^. Figure [l] represents two instances ip and tp of the function 
^ for two different shapes of the room. Let R{S^) denote the group of orientation-preserving rigid motions 
of C M^. We observe that a homeomorphism f : exists, such that ip = o f and / ^ R{S^). It 

follows that djjomeo(s^){^T'4') = 0, so that classical Persistent Homology cannot give positive lower bounds 
for (ifl(5'i)(iy9, -0), while we will see that dii(s^){,p,'4') > 0- 



a 




Room A Room B 

Figure 1: Two rooms and the respective functions ip, i/i, representing minus the distance between the center and the walls. 
is identified with the interval [0, 2-k\. 

Fortunately, we can adapt Persistent Homology in order to obtain a theory that can give a positive lower 
bound for da, in the previous example (and in many similar cases). We are going to describe this idea in 
the next section. 

2. Adapting Persistent Homology to the group G 

Shape comparison is commonly based on comparing properties (usually described by M'^-valued functions) 
with respect to the action of a transformation group. Let us interpret these concepts in a homological setting. 
Before proceeding, let us fix a chain complex (G, d) over a field K (so that each group of n-chains G„ is a 
vector space). We consider the partial order < on E*^ defined by setting [ui, . . . ,Uk) di {vi, ■ ■ ■ ,Vk) if and 
only if Uj < Vj for every j € {1, . . . , k}. 

Definition 2.1. Assume a function (p = (1^1, . . . , pk) ■ Un G„ K*^ U {—00, . . . , —00) is given, such that 
i) (p takes the null chain £ C„ to the /c-tuple (—00, . . . , —00), for every n G Z; 
a) (p{dc) ^ (p{c) for every c e 1J„ G„; 
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in) (p{Xc) = (p{c) for every c € 1J„ C„, A G K, A 7^ 0; 

if) ipj{ci + C2) < max ((^^(ci), (pj{c2)) for every ci, C2 G C„ with n e Z, and every j e {1, . . . , fc}. 
We shall say that <^ is a filtering function on the chain complex (C, d). 

Definition 2.2. Let us assume that a group G is given, such that G acts linearly on each C„ and its action 
commutes with d, i.e., d o g = g o d for every g € G (in particular, every 5 e G is a chain isomorphism 
from C to C). The chain complex {G,d) will be said a G-chain complex. We shall call the group Hn{G) := 
ker9„/im 9„+i the n-i/i homology group associated with the G-chain complex {G,d). 

Now, let us assume that (C, d) is a G-chain complex, endowed with a filtering function (^. For every 
M G R*^ we can consider the chain subconiplex G*^-" of G defined by setting G**-" := {c e G^ : (p(c) ^ u} and 



restricting 9 to G"^^". G*'^" is a subcomplex of G because of the properties in Definition |2.l| (in particular, 
9(Gf|") C G^f^^"). We observe that G^^"" will not be a G-chain complex, since g(G^=^") ^ G;f=^", in 
general. For the sake of simplicity, we will use the symbol d in place of d^(j^^u. 

Definition 2.3. The chain complex (G'^-",(3) will be called the chain subcomplex of {G,d) associated with 
the value u G M.^ , with respect to the filtering function (p. 

We refer to ^5j for the definition of chain subcomplex. 

Now we can define the concept of the n-th persistent homology group of (G, 9), with respect to 

Definition 2.4. If m = (ui, . . . , Uk),v — {vi, . . . , Vk) £ and u < v (i.e., Uj < vj for every index j), we can 
consider the inclusion i of the chain complex G*^-" into the chain complex C^-^ . Such an inclusion induces 
a homomorphism i* : {G'^^'') H„ (G**^"). We shah call the group PH^{u,v) := i* (ff„ (G*^^")) the 
n-th persistent homology group of the G-chain complex G, computed at the point {u,v) with respect to the 
filtering function (p. The rank pf^iu^v) of this group will be said the n-th persistent Betti number function 
(PBNF) of the G-chain complex G , computed at the point (u, v) with respect to the filtering function (p. 

The key property of PH^ is the invariance expressed by the following result. 

Theorem 2.5. If g ^ G and u,v ^ M*"' with u ^ v, the groups PHf^°^{u,v) and PH!^{u,v) are isomorphic. 

Proof. We define a map F : PH^°s{u,v) PH!^{u,v) in the following way. Let us consider an element 
z e PHf^°9{u^ v) := i* {Hn (G'^°9^")) . By definition, a cycle c e G^°f=^" exists, such that z is the equivalence 
class [c]y of c in ff„ (G'^"^-"). We observe that g{c) G G^"-" and the equivalence class [^(c)]^ of g{c) in 
Hn {Gf^") belongs to PH^{u,v) := i* (ff„ (G^^")). We set F{z) = [g(c)]„. 

If c' e G,f 9=^" is another cycle such that z ^ [c']„ e i7„ (G*'°9^"), then a chain 7 e Gf^i-" exists, such 
that c' — c = d-f. We observe that g{j) e G^^^. The inequality (f{d{g{'y))) ^ tp{g{j)) (see Definition 2.1 1 
implies that 9(3(7)) G C*^-"- As a consequence, [g{c')]y = [g{c-\-dj)]y = [g{c)-\-g{dj)]v = [51(c) -H9(5(7))]7^ 
[^(c)]^ + [d{g{j))]v = [5(c)],,. These equalities follow from the linearity of g and the equality d o g = g o d 
in Definition |2.2| This proves that F is well defined. 

Let zi = [ci]y,Z2 = [c2]v e PHf^°s{u,v), with ci,C2 e G^^^". We observe that (7(ci),5(c2) S G,f=^". 
From the linearity of g, it follows that g{\ici -\- A2C2) = Ai(7(ci) + \2g{c2) € G^-"^ , for every Ai, A2 G K. 
Hence, we have that F(AiZi + A22:2) = F(Ai[ci]„ + A2[c2]t,) = i^([AiCi + A2C2]i,) = [g(AiCi + A2C2)]t, = 
Ai[5(ci)]i, + A2[5(c2)]„ = \iF{[ci]y) + A2F([c2]^,) = Aii^(zi) + A2F(z2). Therefore, F is linear. 

Furthermore, if F{zi) = F{z2) then [(7(01)]^ = [9{c2)\v, so that a chain 7 G G^^i exists, such that 
]{ci — C2) = ff(ci) — 5(22) — 97. Moreov er, 9 ~"^( 7) € G^^^-". It follows that Ci — C2 = g^^{d^) 



d (ff"^(7)) G G,f°f=^'', because of Definitions 2.1 and 2.2 As a consequence, [ci]^ = [c2]t). This proves that 
F is injective. 

Finally, F is surjective. In order to prove this, we observe that if w G PH^{u,v) := i* (if„ (G*^-")) 
with the homomorphism i* : H,, (G'^=^") H„ (G'^^^") induced by the inclusion i : G'^=^" G'f^'", 
then a chain c G G,f=^" exists such that w = [c]„ G iJ„ {G^^'"). We have that ^"^(c) G G,?°f=^" and 
J^([5"'(c)].) = [c]v^w. 

Therefore F : PH:^°^{u,v) PHI^{u,v) is an isomorphism. □ 



3 



The previous theorem justifies the name G-invariant Persistent Homology^ showing that the PBNFs of 
a G-chain complex do not change if we substitute the filtering function (f with the function (fog, for g (z G. 

3. Stability of the PBNFs with respect to da 

Let X and {S{X), d) be a triangulable space and its singular chain complex over a field K, respectively. 

Assume that a subgroup G of the group Homeo{X) of all homcomorphisms f : X ^ X and a continuous 
function — ((^i, . . . , Lpk) : ^ — > K*^ are chosen. For every u £ R'^, let us set X*^-" -.^ {x E X : (p{x) ^ u}. 
Let us consider the action of G on S{X) defined by setting 5(17) := g o a for every g (z G and every singular 
simplex a in X, and extending this action linearly on S{X). We recall that, by definition, every singular 
n-simplex in X is a continuous function from the standard n-simplex into X. 

Now, assume that a G-chain subcomplex (G,9) of the singular chain complex {S{X),d) is given. We 
observe that, for every topological subspace X of X, (Gn S{X), d) is a chain complex over the field K. The 
symbol G n S{X) denotes the chain complex G' where G^ is the vector space of the singular n-chains in X 
that belong to G„. 

In order to avoid "wild" chain complexes, we also make this assumption: 

(*) If X' and X" arc two closed subsets of X with X' C int{X"), then a topological subspace X oi X exists 
such that X' C X C X" and the homology group i?„(G n S{X)) is finitely generatedj^ 

Let us consider the set {cr"}jeJ of ^lU (distinct) singular n-simplexes in X. Obviously, ii X is not a finite 
topological space, J will be an infinite (usually uncountable) set. Then we can endow the chain complex G 
with a filtering function (p in the following way. If c equals the null chain in G„, we set if{c) := (—00, . . . , —00). 
If c is a non-null singular n-chain, we can write c = '^'^'^jt ^ " with a'' € K, a*" 7^ for every index r, 

and jr' 7^ jr" for r' ^ r" . In this case we set (p{c) = (ui, . . . ,Uk) € M*^, with each Ui equal to the maximum 
of ipi on the union of the images of the singular simplexes (j"^^ , • ■ ■ , c"^ . In other words, <p(c) is the smallest 
vector u such that the corresponding sublevel set X'-'^-" contains the image of each singular simplex cr"^ 
involved in the representation of c. We observe that this representation is unique up to permutations of its 
summands, so that (p is well defined. Furthermore, the properties in Definition |2.1| are fulfilled. We shall 
say that the function (p is induced by ip. 

An elementary introduction to singular homology can be found in |14j . 

The next result has a key role in the rest of this paper. 

Proposition 3.1. For every n G Z the n-th persistent Betti number function p^{u, v) of the G-chain complex 
(G, d), endowed with the filtering function (p, is finite at each point (u, v) in its domain. 

Proof. Since u ^ v and ip is continuous, we have that the set X^-"^ is closed and contained in the interior of 
the closed set X"^-"". Property (*) implies that a topological subspace X of X exists such that X''^-" C X C 

X^^" and Hn{C n S{X)) is finitely generated. The inclusions G n S'(X'^^«) A G n S{X) 4 G n SiX^^"") 

induce the homomorphisms7J„(Gn5(X^=^")) ^ Hn{Cr\S{X)) ^ Hn{Cr\S{X'^^'")). Since_dimim (j*oi*) < 
dimimj* < dim i?„(G n S{X)) < +00, we obtain that also PH^{u,v) := j* o i* (iJ„ [G f] S {X'f^''))) is 
finitely generated. □ 

From now on, in order to avoid technicalities that are not relevant in this paper, we shall consider two 
PBNFs equivalent if they differ in a subset of their domain that has a vanishing measure. 

A standard way of comparing two classical persistent Betti number functions is the matching distance 
Dmatch, a.k.a. bottleneck distance (cf. [9l[5]). It can be applied without any modification to the case of the 
persistent Betti number functions of the G-chain complex G. An important consequence of the finitcness 



^We wish to avoid chain complexes like the one where the 0-chains are all the usual singular 0-chains and the only 1-chain is 
the trivial one. In this case the homology group Hq{C) would not be finitely generated, in general. This means that property 
(*) would not hold for X' = X" = X. 



4 



of these functions is the fohowing theorem, showing that the matching distance between persistent Betti 
number functions of the G-chain complex C is a lower bound for the natural pseudo-distance da- In other 
words, a small change of the filtering function with respect to da produces just a small change of the 
corresponding persistent Betti number function with respect to D^atch- This property allows the use of 
PBNFs in real applications, where the presence of noise is unavoidable. 

Theorem 3.2. For every n G Z, let us consider the n-th persistent Betti number functions , pfl of the G- 
chain complex {C, d), endowed with the filtering functions (p and 'i/j induced by (p : X Mf' and ij) : X M'^' , 
respectively. Then Dmatchipti pt) < <^g('^,'0)- 

Proof. We can proceed by mimicking the proof of stability for ordinary persistent Betti number functions 
(cf. [5j). This is possible because that proof depends only on properties of PBNFs that are shared by both 
classical persistent Betti number functions and persistent Betti number functions of a G-chain complex 



endowed with a filtering function, once we have proven that the PBNFs are finite (Proposition 3.1). It is 
sufficient to substitute the group Homeo{X) with the group G C Homeo{X), and the homology groups of 
each sublevel set X'-''-" with the homology groups of the G-chain complex C n S'(X'^-"). □ 

4. Applications 

In this section we illustrate how G-invariant persistent homology can be used to discriminate between 



the rooms described in Example 1.2 showing that no rotation of changes the function tp into ip. 

In order to manage this problem we can consider the i?(S'^ )-chain complex C whose 7i-chains are the 
singular n-chains c S Sn{S^) for which the following property holds: 

(P) If a singular simplex cr" appears in the representation of c with respect to the basis {c"} of Sn{S^), 
then the antipodal simplex s o cr" appears in that representation with the same multiplicity of cr", 
where s is the antipodal map s : ^ . 

In other words, we accept only chains that can be written in the form X^^Li o,^ ('^jt + ° "^jv) ■ E^^^'Y rotation 
p £ R{S^) commutes with the antipodal map s and is a chain isomorphism from C to G. Moreover, it is 



easy to verify that the properties in Definition 2.2 are fulfilled, for G = R{S^) and C — C. The chains in G 
will be called symmetric chains. 

We can prove that property (*) holds for the i?(5'^ )-chain complex that we have defined. Let X' and X" 
be two closed subsets of with X' C int{X"). Let us set X equal to the e-dilatiorrl of X' in S*^, choosing 

£ > so smah that the X C int{X"). We observe that the set X{^s{X) is open and s n s{X)^ = X{^s{X). 

Moreover, Xr\s{X) is the union of a finite family J- = {ai} of pairwise disjoint open arcs, having the property 
that iiai € F then also s{ai) S F (possibly, F = {S^}). Now, let us consider the topological quotient space 
Q obtained by taking all unordered pairs of antipodal points in X H s{X). We have that Q is homeomorphic 
to the union of a finite family F' of pairwise disjoint open arcs of (possibly, F' = {"S*^}), and hence the 

ri-th homology group Hn{Q) is finitely generated. A chain isomorphism F from Gn 5 [xc^s{X)^ to S{Q) 

exists, taking each chain cr + s o cr to the chain given by the singular simplex {a, s o cr} in Q. F induces an 

isomorphism from Hn (c H S (^X n s(A))) to HniQ). Therefore also Hn{Cf^s{xr\ s(A)^^ is finitely 

generated. Property (*) follows by observing that G n (a n s(A)) = G n S{X). 

Referring to Example |1.2[ we observe that the matching distance between the 0-th persistent Betti 
number functions of the i?(5'^ )-chain complex G with respect to the filtering functions (p and is positive. 
Hence, Theorem |3.2| gives a non-trivial lower bound for djf(gi'f{ip,ip), while the matching distance between 
the corresponding classical persistent Betti number functions vanishes. The previous claim becomes clear 



^The e-dilation of a subset y of a metric space M is tlie set of points of M that have a distance strictly less than e from Y. 
On C we consider the metric induced by the Euchdean metric in R^. 
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if we consider the birth of the first homology class in the homology groups Hq (C"^-*) and Hq ^C"^-*^, 

respectively, when the parameter t increases. While the group Hq (C"^-*) becomes non-trivial when t 

reaches the value to = imn ip ~ min?/;, the group Hq becomes non-trivial when t reaches a value 

i > mm ip = mmtp. This is due to the fact that the sublevel set {a; € S*^ : ip{x) < to} contains two pairs 
of antipodal points, while the sublevel set {x G : ipi^) ^ ^o} contains no pair of antipodal points (see 



Figure^. By applying Theorem 3.2 it follows that dfl(5i)((/3, tp)>t — Iq. 



The interested reader can find the 0-th persistent Betti number functions p*^ and p'f^ of the i?(5'^)-chain 
complex C in Figure [sj 



These points 
are antipoda 



These points 
are antipodal 




These points These points 
are antipodal are antipodal 




■ / 1 \ ' - n/ 



Figure 2: The sublevel sets of the filtering functions ip^ i/i cited in Example |1.2[ respectively for the levels to and t. 




Figure 3: The 0-th persistent Betti number functions pj^ and of the H{5^)-chain complex 0, corresponding to the filtering 
functions i/p, -0 cited in Example |1.2| In each part of the domain, the value taken by the PBNF is displayed. Observe that in 
both figures a small triangle is present, at which the persistent Betti number function takes the value 2. 



Remark 4.1. As an alternative approach to the problem of comparing two filtering functions ip,ip : X R, 
the reader could think of using the well known concept of Equivariant Homology (cf. |16j). In other words, 
in the case that G acts freely on X, one could think of considering the topological quotient space X/G, 
endowed with the filtering functions ip, ip that take each orbit uj of the group G to the maximum of p> and "0 
on Lu, respectively. We observe that this approach would not be of help in the case illustrated in Example |1. 2 [ 
since the quotient of / R{S^) is just a singleton. As a consequence, if we considered two filtering functions 
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(p,ip ■ ^ R with max ip — maxi/), the persistent homology of the induced functions (p, ip ■ /R{S^) — > M 
would be the same. 
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